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Abstract 

The vertical component (Z) of the geomagnetic field observed by ground-based observatories of 
the INTERMAGNET network has been used to analyze the effects of the movement of electri- 
cally conducting sea water through the geomagnetic field due to a propagation of a tsumani. The 
purpose of this work is to study the geomagnetic variations induced by the tsunamis occurred at 
26 December, 2004, 27 February, 2010 and 11 March, 2011. For each case study, we selected 
four magnetic stations belonging to the INTERMAGNET programme that were influenced or 
more direct affected by the tsumani. To detect these disturbances in the geomagnetic data, the 
discrete wavelet technique have been used in four levels of decomposition. We were able to 
detect the localized behavior of the geomagnetic variations induced by the movement of elec- 
trically conducting sea-water through the geomagnetic field, i. e., the identification of transients 
related to the tsunamis. As well, using the minutely magnetogram data, it was able to localize 
the initial phase and time of the tsunami maximum height. The first interpretation of the results 
suggests that discrete wavelet transform can be used to characterize the tsumani effects on the 
geomagnetic field, but need further study. 
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1. Introduction 

In the treatise "De Magnete" as early as 1600, the english researcher Gilbert showed the 
predominately dipolar character of the terrestrial magnetic field. Next in the early nineteen cen- 
tury, other english researcher Gauss improved the magnetic field observation techniques and 
introduced the spherical harmonic method for geomagnetic field analysis. Based on a called 
"Spherical Harmonic Analysis", a mathematical method was developed to separate the external 
and the i nternal contributi ons to the surface field by a unique global analysis of the Earth's main 
field (see lCampbel][ll997l and references therein). 



The geomagnetic field is described as a complicated function of space and time. Ground 
based magnetic measurements show a repetitive diurnal variation on geomagnetically quiet days. 
But there is a great variety of irregular variations that occur from time to time, the "disturbance 
fields". Periods o f great disturbance are called, by analogy with the weather, "magnetic storms" 



(lParkinsonLll983l) . 

Some e vidence of the influenc e of oceanic tides on the magnetic daily variation has been 
obtained by Larsen and Cox ( 19661) . They found small semidiurnal variations of the Z component 



at a coastal site (Cambria, California) and at two island stations (Honolulu and San Miguel) that 
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could not be explained by the atmospheric tidal theory. They suggested that these variations 
must be due predominantly to oceanic tides. It is important to mention here that the conductivity 
of the ocean does not vary significantly with time, unlike the ionospheric conductivity. As a 
consequence, the seasonal variation of the oc eanic contribution is expected to be smaller than the 



ion ospheric con t ributio n (ICueto et al.U 2003). 



Manoi et al. (2011) observed the geomagnetic contributions due to the moderate tsumani 



in the Pacific ocean provoked by the Chilean earthquake of 8.8 in magnitude at three different 
magnetic stations (PPT, HUA and IPM). Their observations presented a variation of 1 nT in 
the vertical component of the magnetic field (Z) during the time of the tsumani effects. Move- 
ment of electrically conducting sea-water through the geomagnetic field generates an electromo- 
tive force that induces electric fields, currents and secondary magnetic fields. In other words, 



28 tsunamis 



can produce p erturbations in the Earth's magnetic field by electro-magnetic induction 



29 (see lManoi et all 1201 ll and references therein) 



In this work, we focused on the survey of geomagnetic variations induced by some selected 
tsunamis: one of 26 December 2004, other of 27 February 2010 and another of 1 1 March, 201 1. 
As t he wavelet transforms h ad turn out to be a very useful tool in atmospheric signal analysis 



(see iDorningues et all I2005L and references therein), the discrete wavelet technique have been 



selected and used with four levels of decomposition in order to detect the small singularities 
in the geomagnetic data. For a good revision on the application o f dis crete wavelet on Space 
Weath er Research, Domingue s et al. (2005), iMendes et al. (2005) and Mendes da Costa et al 



d201 lb can be consulted. Thus, this work aims to verify the use of the wavelet technique as a way 
to identify the effects related to tsunamis on the geomagnetic field components, particularly the 
Z-component. The contents are composed by: in Section|2] the presentation of the geomagnetic 
disturbance mechanism related to the tsunami effects; in Section [3] applied the methodology; 
in Section H] the data; in Section [5] the results achieved by the analysis; and in Section |6j the 
conclusions of our work. 



2. Geomagnetic disturbance mechanism 

Faraday originally predicted that ocean waves could induce electrical fields, which are also 
called by "motional induction". The induction mechanism process consists of the flow of the 
ocean water, as an electrically conducting fluid, through the Earth's main magnetic field. The 
salts dissolved in the sea water are made of ions that can be deflected by the Lorentz force which 



48 acts in a perpendicular direction both to the velocity and magnetic field vectors ( Tyl er et al 



2003l:lManoi et al.U 2006). As consequence of the motion of sea water with velocity V across the 



geomagnetic field B, a conduction current J is produced expressed by: 

1 = cr(£ + V XB), (1) 

— > — > — > 
where <x is the electrical conductivity, E is the stationary electric field and V x B is the induced 

electric field. Discarding any displacement current, the Ampere's law is used to determine an 
— > 

induced magnetic field b given by: 

Vxl?=ju ~7, (2) 
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where fi a = 47rl0 7 henry /m. It is known that electric and magnetic fields are related by 
Faraday's law given by: 



■3 -» db 
V x E = — . 

dt 



(3) 



Applying the same mathematical development used by Podnev (1975), the velocity field 
present in the problem can be described as a curl of a vector stream function *F (the stream 
understood as the flow of the ocean water), i.e.: 



V = V x W. 



(4) 



In analogy, the induced magnetic field can be expressed by a vector potential A, given by: 



V X A. 



(5) 



Thus the relation between the vector stream function *P and the geomagnetic field B as a function 
— > 

of the vector potential A can be established as: 



V 2 A -n ct — = -Ho or V X (T X B). 
at 



(6) 



In a short description, the induced magnetic field generated by the ocean can be classified by 
two components: toroidal and poloidal. The toroidal component is generated by the electrical 
currents closing in the vertical plane and can reach up to 100 nT but is confined to the ocean and 
the upper crust. The poloidal component is much weaker, between 1 and 10 nT, and arises from 
the electrical currents closing in the horizo ntal plane, however it can reaches out side of the ocean 
to remote lands and satellite locations (see Tyler et al. , 2003 : Manoi et al. , 20061 and references 
therein). 

The sea water is subject to two kinds of geomagnetic fluctuations: peri odic (lunar and solar 
variati ons) and transient disturbances (geomagnetic storms), as discussed bv lLonguet -Higgins et al 
d 1954b . Also, the ocean flow can be occasionally disturbed by tsunamis that can increase the 
ocean wave speed for few centimeters per second and that affects the entire water column. Con- 
sequently, as demonstrated before, to increase the ocean flow will also increase the magnetic 
induced field. 



Considering quiet conditions. iTvler et al.l y999) discussed that the major difficulty to deter- 
mine the magnetic ocean generated signals are the weak values compared to the signals from 
other sources. This is the reason to take advantage of some analysis features of the wavelet 
technique. 



3. Discrete Wavelet Transform (DWT) methodology 



The wavelet transforms were better and bro adly formalized thanks to mathematicians, physi- 
cist, and engineers efforts (e.g., iMorlen, 1 1 983b . In the domain of geophysics applications, the 
main c haracteristic of the wavelet technique i s the introduction of time-frequency decomposition 



(e.g.. lDomingues et all 12005: Men des et al.L 2005; Me ndes da Costa et 



requency decomposition 
ail 1201 ll) . In the 1990s 



several important i de as and applications concerning wavelet were developed (e.g., Daubechies, 
119921: IChuiL fl992a1lrl 1 1 994h . 
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The wavelet analysis could show that the larger amplit udes of the wavelet coefficients are lo- 
cally associated with abrupt signal ch anges. In recent work (iDomingues et all 2005; M endes et al 



2005 : Mendes da Costa et al. , 201 ll) . a method for the detection of the transition region and the 



exact location of this discontinuities due to geomagnetic storms was implemented. 

As a basis for this work, this mathematical method is presented in an abbreviated way. The 
multiresolution analysis (AMR) is a mathematical tool used to build up wavelet functions (e.g., 
Dau bechiesl 1 19921: iMallatl Il99lh . Then it is possible to build up wavelet functions using a pair 
of vectors { W, <pi}, in such a way that there are a sequence of embedded approximating spaces 
and the functions W c V /+1 and (p J k form a Riesz basis for V J , being V J = span{(/> J k (t)}. In this 
technique, a mother-wave let is genera ted from a scaling function. It obeys the following scale 



relation (more details in Mallat, 199 lh : 



(f,(t) = 2 V h(k)<t>(2t - k) (Scale relation) 



(7) 



keZ 



and h(k) is the low pass-filter, also called scale filter coefficients. The family of scale relations 
are represented by: 

(8) 



cp ] k (t) = 2 j/2 (p(2 J 't - k), j,keZ 

99 They are called scale functions and in frequency domain are given by: 

100 where H(£) = £ e~ lk ^ is a low pass filter associate to <p. The following relation holds: 



ke h(k) 



V j+1 = V j e W j 



(9) 



(10) 



The spans W J have the difference of information between V and V 7+1 . If the \p function form a 
Riesz base of W J , it is called wavelet function: 



= 2«hK2h - k) m = 2 Yj g(W(2t - k) 



(11) 



keZ 



where g(k) = (-1) • h(-k + 1) is a high pass filter. 

The wavelet and scale functions satisfy the orthogonality condition, (<p J k , ifrj) = 0, {tyl, ij/J 1 ) = 
dj, m 5k,i- The AMR tool is useful to study the function in L 2 (!R). The difference of information 
between V-i and V /+I is given by i. e., 



(W +1 -W)f(f) = Qjf(t), 
where the projections in W and W-* are, respectively: 

n j f(t) = 2</, <f>i)4(t) Q j f(t) = ifi)ifi(f), 

k k 

We obtain in multi-level jo < j the coefficients expression: 



(12) 



(13) 



(14) 
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We can obtain with a change of base \(f>l +l } 



U 



. U {iffi}, the equation: 



In the DWT the equation (fT4t is manipulated jointly with to the scale relations: 
c[ = V2 J] A(m - 2^)c4 +I = V5 £ g(m - 2k)c^ 



(15) 



(16) 
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In this work we use the Daubechies (dbl) wavelet function of order 1 (Haar wavelet). In this 
case the coefficients h = [-^, -^] and consequently, g - [-^, — -^]. 

The wavelet transform can be used in the analysis of non-stationary signals to rescue infor- 
mation on the frequency or scale variations of those signals and to detect the localization of theirs 
structures in time and/or in space. Due to the double localization property of the wavelet func- 
tion, the wavelet transform is said to be of local type in time-frequency , with time and frequency 
resolutions inversely proportional. As discussed by iDomingues et al.l (120051) . one consequence 
of the properties of the wavelet analysis is the possibility of showing that the wavelet coefficients 
amplitude is associated wi th abrupt signal variations or "details" of higher frequency (see also 
Meven, 1 1 990t IChuil 1 1 992bl) . Selecting a wavelet function which closely matches the signal to be 
processed is of utmost importance in wavelet applications. A tutor ial on the main properties and 
applications of the w avelet analysis can be found, for instance, in IDomingues et al.l (12005b and 
Mendes et al.l d2005h . 

A hard thresholding process has been applied to the wavelet coefficients for the first three 
decompositio n level s to identify the storm time in the magnetograms, as has been discussed by 
Mendes et al.l (120051) . Thresholding concerns the process of setting to zero certain coefficients 
in an effort to highlight significant information, in this c ase the shock-like transient phenomena 
that a ppear in the main phase of the magnetic storms (see Domingues et al. , 2005 ; Mendes et all 
2005b . It was found that could be used as threshold sets for the decomposition levels (d- 7 ) 2 , 
because accor ding to Meyer theorem, the multilevel threshold sets may have a multiplicity of 2 
dMeveilll99d) . 

In this analysis, we choose to use the Haar wavelet and the sampling rate of 1 min, as con- 
sequence of the time resolution of the magnetic data, with the pseudo-periods of the first four 
levels 2, 4, 8 and 16 min. The Haar wavelet is considered the most simple orthogonal analyzing 
wavelet function. However the discrete wavelet transform (DWT) using Haar wavelet can detect 
abrupt variations, i. <?., one localization feature in the physical space. The Haar wavelet has a 
compact support in physical space and a large support in Fourier space. In other words, Haar 
wavelet has a better localization on ti me than on frequen cy, what is expected by the Heisenberg's 
uncertainty principle (as discussed in Daubechies, 1992). 



4. Magnetic Data 

In this section, we first describe the data used to study the geomagnetic variations due to the 
tsunami-generated magnetic fields. The three tsunami's events are presented. For each event, 
we have chosen four ground magnetic measurements. We have also selected the tide-gauge 
measurements at or nearby the chosen magnetic stations, only for guiding purposes. 
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We selected magnetic stations belonging to the INTERMAGNET programme (http://www 
intermagnet.org) that were influenced or more direct affected by the tsumanis. By international 
agreement, there are two systems t hat can represent the Earth's magnetic field: the XYZ and the 
HDZ system (see lCampbell , 1 19971 and references therein). The X, Y and Z stand for northward, 
eastward and vertical into the Earth directions, the H, D and Z stand for horizontal component, 
declination (angular direction of the horizontal component) and vertical (into the Earth). The H- 
component is more affected by the solar-magnetospheric interactions, consequently, also the X- 
and Y-component. These variations, specially those associated to the ring current, have a major 
contribution on the magnetograms at stations located at low and mid-latitude regions. Due to 
the Z-component be less affected than the H-component, we decided to use the Z-component to 
detect the geomagnet ic variations i nduce d by the tsunami. Using a different technique in a work 
with similar purpose, Ma noi et alJ yOl 1) also used the Z-component to observe the variations in 
the Z-component provoked by the Chilean tsumani, 27 February, 2010. 

The intensity of the geomagnetic disturbance in each day is described by indices. There are 
different indices that can be used depending on the character and the latitude influences in focus. 
Kp, AE and Dst and their derivations are the most used geomagnetic indices. The Kp index is 
a number from to 9 obtained as the mean valu e of the distu rbance levels within 3-h interval 
observed at 13 subauroral magnetic stations (see lBartelsLll957h . The minutely AE index (some- 
times 2.5 minute interval) is obtained by the H-component measured from the magnetic stations 
located at auroral zones. The index most used in low and mid-latitudes is the Dst index. It rep- 
resents the variations of the H-component due to changes of the ring current. The characteristic 
signature of a magnetic storm is a depression in the horizontal component of the Earth's magnetic 
field H and this depression is shown by the Dst index. 

In this work, we analyzed the SYM-H index to identify periods of disturbed geomagnetic 
field conditions. The SYM-H is essentially the same as the traditional hourly Dst index. The 
main characteristic of the 1 minute time resolution SYM-H index is that the solar wind dynamic 
pressure variation are more clearly seen than indices with lower time resolution. Its calculation 
is based on magnetic data provided by 1 1 stations of low and medium latitude. Only six of 
the stations are used for its calculation of each month, some stations can be replaced by others 
depending on the data conditions. 

We focused on the survey of geomagnetic variations induced by the tsunamis occurred at 26 
December, 2004, 27 February, 2010 and 11 March, 2011. Related to the magnetic disturbance 
level, these events of tsunami occurrence happened during quiet time (27 February, 2010) and 
disturbed geomagnetic conditions (weak - 26 December, 2004, and moderate storms, 1 1 March, 
2011). 

On 26 December 2004 the Sumatra- Andaman earthquake occurred with 9.1M W (M w means 
moment magnitude of the Richter scale) at 00 : 59 UT (Universal Time) . The epicenter coor- 
dinates are Lat. 3.4° and Long. 95.7° at 30 km depth. iTsuii et alJ (2006) carried out tsunami 
height measurements on the Andaman seacoast. They found the maximum tsunami height was 
19.6 m at Ban Thung Dap (Lat. 9.03° and Long. 98.26°) and the second highest was 15.8 m at 
Ban Nam Kim (Lat.8.86° and Long. 98. 28°). The tsunami not only propag ated throughtout the 
Indian ocean but also into the Pacific and Atlantic oceans. IMurtv et alJ ([2005) discussed through 
a simple analytical model that the tsunami flux from the Indian ocean into the Pacific ocean is 
greater than the flux into the Atlantic ocean. This leakage of tsunami energy into the Pacific and 
Atlantic oceans occurs through the gaps between Australia and Antarctica and between Africa 
and Antarctica, respectively. The maximum tsunami amplitudes were o f 0.65 m on the Pacific 
coast and 0.3 m on the Atlantic coast as discussed bv lMurtv et al.1 Q2005L and references within). 
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192 To study the geomagnetic variations due to the Sumatra-Andaman tsunami, the magnetic 

193 stations considered were: Eyrewell (EYR), Huancayo (HUA), Learmonth (LRM) and Papeete 

194 (PPT), with the geographic and geomagnetic coordinates presented in TableQ] Figure [4] display 

195 the magnetic stations distribution and their IAGA codes. The selected tide-gauge measurements, 

196 only for guiding purposes, are located at Cocos Island (Lat. -12.12°, Long. 96.90°), Jackson 
,97 Bay (Lat. -43.97°, Long. 168.62°), Nuku Hiva (Lat. -08.82°, Long. 140.20°) and Callao (Lat. 
,98 -12.08°, Long. -77.13°). 



Table 1 : INTERMAGNET network of geomagnetic stations for the study of the Thai tsunami. 

Station Geografic coord. Geomagnetic coord. 



Lat.(°) Long.Q Lat.(°) Long.(°) 

EYR -43.41 172.35 -46.92 -106.22 

HUA -12.04 -75.32 -1.73 -3.44 

LRM -22.22 114.10 -32.28 -173.53 

PPT -17.75 -149.57 -15.20 -74.83 

Source: http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html (2005) 




Figure 1 : Geographical localization of the stations used in the study of the Sumatra- Andaman tsunami 2004 

199 Near coast of the central Chile, on the 27 th of February, 2010 at 06 : 34 UT occurred 

200 an earthquake with magnitude 8.8 M w . The epicenter was loca ted on Lat. -36.1° and Long. 



213 



-72.6° at 55 km depth. As reported bv lPararas -Caravannisl (2010). shortly after the earthquake, 



202 tsunami waves hit the coastal area of the Central Chile. The tsunami overtook the coastal cities 

203 as Talcahuano, Coquimbo, Antofasta and Caldera, as well as the Juan Fernandez Islands. The 

204 NOAA Pacific Warning Center released a bulletin number 018 and a tsunami warning was issued 

205 at 00 : 12 UT on 28 February, 2010 for Chile, Peru, Ecuador, Colombia, Antarctica, Panama, 

206 Costa Rica, Nicaragua, Pitcairn, Honduras, El Salvador, Guatemala, French Polynesia, Mexico, 

207 Cook Islands, Kiribati, Kermadec Island, Niue, New Zealand, Tonga, American Samoa, Samoa, 

208 Jarvis Island, Wallis-Futuna, Tokelau, Fiji, Australia, Hawaii, Palmyra Island, Johnston Island, 

209 Marshall Island, Midway Island, Wake Island, Tuvalu, Vanuatu, Howland-Baker, New Cale- 

210 donea, Solomon Island, Nauru, Kosrae, Papaua New Guinea, Pohnpei, Chuuk, Marcus Island, 

211 Indonesia, North Marianas, Guam, Yap, Belau, Philippines and Chinese Taipei. 

212 The magnetic stations considered for this event were: Huancayo (HUA), Easter Island (IPM), 
Papeete (PPT) and Teoloyucan (TEO), with the geographic and geomagnetic coordinates pre- 



sented in Table|2] Figure[4]display the magnetic stations dis tribution and their IAGA codes. We 
selected the same three stations used bv lManoi et al. ( 201 ll) to study the geomagnetic contribu- 
tions due to the Chilean tsumani and included the station of TEO located in Mexico. In their 
study only the IPM station showed a periodic variation of 1 nT in the vertical component (Z) 
caused by the tsunami started at 1 1 : 35 UT, and the other two stations did not show concurrent 
variations. 

The selected tide-gauge measurements are located at Callao La Punta (Lat. -12.01° and 
Long. -77.17°), Papeete (Lat. -17.75° and Long. -149.57°) and Acapulco (Lat. 16.83° and 
Long. -99.92°). 



Table 2: INTERMAGNET network of geomagnetic stations for the study of the Chilean tsunami. 



Station 


Geografic coord. 


Geomaj 


;netic coord. 


Lat.(°) 


Long.(°) 


Lat.(°) 


Long.(°) 


HUA 


-12.04 


-75.32 


-1.99 


-3.03 


IPM 


-27.90 


-109.25 


-19.63 


-34.47 


PPT 


-17.75 


-149.57 


-15.20 


-74.49 


TEO 


19.75 


-99.19 


28.45 


-29.07 



Source: http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html (2010) 




Figure 2: Geographical localization of the stations used in the study of the Chilean tsunami, 2010. 

223 On 1 1 March, 201 1, an earthquake occurred in the Japanese coast at 05 : 46 UT with 8.9 M w . 

224 The epicenter located at Lat. 38.3 ° and Long. 142.4°, near to the east coast of Honshu, at 

225 24 km depth. That earthquake induced a "major" tsunami that affected all the Japanese coast and 

226 transversed the Pacific Ocean eastward. 

227 As we have done before , we also choose four stations that belong to the INTERMAGNET 

228 programme that were influenced or more directly affected by the Japanese tsunami of March 11, 

229 201 1 to study the geomagnetic variations. The stations considered in this analysis were: Charters 

230 Towers (CTA), Guam (GUA), Kanoya (KNY) and Memambetsu (MMB), with geographic and 

231 geomagnetic coordinates presented in TablefJ] Figure|4]display the magnetic stations distribution 

232 and their IAGA codes. 
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233 The selected tide-gauge measurements are located at Hanasaki (Lat. 43.28° and Long. 

234 145.57°), Tosashimizu (Lat. 32.78° and Long. 132.92°), Pago Bay (Lat. 13.43° and Long. 

235 144.80°) and Cape Ferguson (Lat.-19.28° and Long. 147.06°). 



Table 3: INTERMAGNET network of geomagnetic stations for the study of the Japanese tsunami. 



Station 


Geografic coord. 


Geoma| 


;netic coord. 


Lat.(°) 


Long.(°) 


Lat.(°) 


Long.(°) 


CTA 


-20.10 


146.30 


-27.59 


-138.58 


GUA 


13.28 


144.45 


05.18 


-144.20 


KNY 


31.42 


130.88 


22.06 


-158.71 


MMB 


35.44 


144.19 


27.09 


-146.95 



Source: http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html (2010) 




Figure 3: Geographical localization of the stations used in the study of the Japanese tsunami 201 1 



236 5. Results and analysis 

237 Initially the wavelet treatment is done on the SYM-H Dataset, which is a processed data. 

238 According to the theory used, expected results are obtained. After that the same treatment is 

239 done directly on the magnetometer dataset for each tsunami event. The graphical results are 

240 presented and discussed validating the identifications of the tsunami occurrences (Dec. 26, 2004, 

241 Feb. 27, 2010, and Mar. 1 1, 201 1). 

242 From theory the wavelet coefficients are supposed to allow identifying the induced magnetic 

243 variations due to tsunamis, because the technique can detect physical discontinuities in the ver- 

244 tical component of the geomagnetic field. When the magnetosphere is under quiet conditions 

245 the behavior of the recorded Z-component should be much smoother than its behavior in the dis- 

246 turbed periods, making easy the identification of the variations only induced by the propagation 

9 



247 of the tsunami. Using a very localized wavelet, we will be able to detect these small geomag- 

248 netic variations. The Haar wavelet present a short compact support that provides a better local 

249 characterization of the signal and it may be the most convenient to represent times series with 

250 abrupt variations or steps. 



252 
253 
254 
255 
256 
257 
258 
259 
260 



5.1. Analysis on SYM-H Dataset 

Figure |4] presents the graphical results of wavelet analysis for sample periods related to 
tsunami occurrences. In each part, at the top there is the SYM-H record and at the bottom 
respectively from the first to the fourth level of the wavelet decomposition. Based on the theory, 
the amplitude of the wavelet coefficients should identify the transient behavior in the dataset. 

In the figure it is shown the discrete wavelet decomposition applied to minutely SYM-H 
index using Daubechies orthogonal wavelet family 1 . From top to bottom each panel shows, the 
SYM-H and the first four levels of wavelet decomposition {(d^) 2 where j = 1,2,3,4) for the day 
corresponding to the tsumani of 2004 due to the Sumatra- Andaman earthquake (Figure |4^) and 
the day after (Figure |4j)). 



j -50- 
M -100- 



(a) 26"' December, 2004 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 
Hours [LIT) 



P 50 



100 — 
S 50- 



1000- 
5 500- 



(b) 27 f/! December, 2004 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 
Hours (UT) 



Figure 4: Graphics of SYM-H with the corresponding wavelet coefficients (d^) 1 for j = 1, 2, 3, 4. The panel corresponds 
to a tsumani occurrence, (a) the arrival day and (b) the day after the Sumatra- Andaman tsunami 2004. 

261 The SYM-H index analysis for the day 26 December, 2004 showed a maximum positive 

262 value of 8 nT during the period of 03 : 26 UT and 03 : 44 UT and again at 14 : 53 UT and 

263 the minimum negative value of -5 nT occurred at 17 : 00 UT (see Figure^). In first level of 

264 decomposition, (d^) 2 (j = 1) showed two maximum of wavelets coefficients at 08 : 10 UT and at 

265 10 : 38 UT. In second level, (d j ) 2 (j = 2) presented the three highest values at 07 : 28 UT, 08 : 20 

266 UT and 11 : 44 UT. Only one maximum of wavelets coefficients {d^) 2 (j = 3,4) are shown in 

267 the third and fourth levels at 07 : 28 UT and 01 : 36 UT, respectively. Geomagnetic activity is 

268 classified by "size" (the amplitude of the variation in the magnetic records) and usually described 

269 by the variation of indices to distinguish between a quiet and an active day (occurrence of storm 

270 or substorm). The characteristic signature of a magnetic storm is a depression in the horizontal 
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component of the Earth's magnetic field H at middle to low latitude ground stations. The index 
most used in low latitudes is the Dst index, similar to the SYM-H index. However SYM-H is a 
minutely index, while Dst is a hourly index. This is the reason that the wavelet treatment were 
done on SYM-H dataset. They represe nt the variations of the H component due to changes of 



the ring curr ent ( Gonzalez et all 1994 ). Following the terminology of Sugiura and Chapman 
discussed by Gonzalez et alJ ( 1994f) . the 26 December, 2004 presented a minimum Dst value of 



-22 nT at 19 : 00 UT which corresponds to a weak geomagnetic storm. As showed on Figure^, 
this day did not presented abrupt variations and consequently, in all the decompositions levels, 
the wavelet coefficients were quite not big enough. 

We also analyzed the SYM-H index of the day after the Sumatra-Andaman tsunami. On the 
27 December, 2004, the SYM-H presented the maximum positive value of 14 nT at 14 : 55 nT 
and the minimum negative value of -13 nT at 21 : 36 UT and at 23 : 44 UT. The Dst in- 
dex showed minimum negative value of -23 nT between 22 : 00 UT and 24 : 00 UT which 



corresponds to a weak geomagnetic storm. iMendes et alJ d2005l) discussed that when the magne 



285 tosphere is under quiet conditions or the analyzed function is smooth, the wavelets coefficients 

286 show very small amplitudes compare to when a geomagnetic storm is under development or 

287 abrupt functions. The first and second level of decomposition showed the maximum of wavelets 

288 coefficients (d j ) 2 (j = 1,2) at 09 : 16 UT and at 09 : 16 UT, respectively. Both the third and 

289 fourth level showed the maximum of wavelets coefficients {d^) 2 (j = 3,4) at 07 : 12 UT (see 

290 Figure^). The wavelet coefficients also were quite not big enough. 

291 As the results are referred to the analysis of a pos-processed magnetic dataset (SYM-H), the 

292 wavelet coefficients were expected to be quite not big enough. The processing of the magnetic 

293 records in order to create a geomagnetic index acts as a filtering process that eliminates discon- 

294 tinuities. So, a result of non identifying tsunami features was already expected. In the next steps 

295 we are going to deal directly with the records (magnetograms) of the selected magnetic stations. 

296 5.2. 26-December-2004 tsunami event 

297 In Figure [5] panels from (a) to (d) correspond to magnetic stations of LRM, EYR, PPT and 

298 HUA, respectively. Each panel, from top to bottom, displays the corresponding tide gauge mea- 
surement (sea level), the SYM-H index (both have been included for comparison purposes), the 
corresponding magnetogram (Z component), and the wavelet signatures established by the square 
wavelet coefficients for the four first decomposition levels, (d-*) 2 for j = 1,2,3,4. The dashed 
rectangle is related to the period of occurrence of the tsunami identified by the methodology 
presented in this work. 

The initial tsunami wave at the Cocos Island atoll (Australia) arrived at 03 : 17 UT and it was 
measured with amplitude of up to 0.33 m. After the first wave arrived, following tsunami waves 
continued for hours. As reported by the Department of Transport - Government of Western Aus- 
tralia (http://www.transport.wa.gov.au/imarine/19383.asp), the initial wave arrived on Exmouth 
(Lat. -21.93°,Long. 114.13°), Australia around 06 : 30 UT. It is expected that the magnetic field 
induced by the arrival of the tsunami will be detected at LRM around 06 : 30 UT and forward. 
Unfortunately, only Cocos Island was the near city to the magnetic station of LRM that we were 

311 able to acquire the tide gauge measurement. 

312 Figure [5^ shows for guiding purposes the sea level for Cocos Island. This was the near- 

313 est station measurement provided by the National Oceanic and Atmospheric Administration 

314 - NOOA (http://wcatwc.arh.noaa.gov/about/tsunamimain.php). We detected a small amplitude 

315 coefficients at 04 : 08 UT, 17 : 00 UT and 20 : 28 UT, in the four first decomposition levels. 
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316 These {d^) 2 {j = 1,2, 3,4) could not been explained by disturbed geomagnetic activity. Also 

317 (d j ) 2 (j = 1, 2, 3) presented small amplitude coefficients at 14 : 38 UT and 22 : 40 UT. 

318 In Jackson Bay, New Zealand, the initial wave arrived at 15 : 10 UT with amplitudes up to 

319 0.65 m. In order to verify the transients due to the tsunami we used the EYR magnetogram data 

320 (see Figure^). For (d- 7 ) 2 (j = 1), the highest wavelet amplitudes were detected at 03 : 58 UT, 

321 05 : 08 UT, 18 : 00 UT, 20 : 36 UT and 21 : 22 UT. We are only detecting the highest amplitudes 

322 that do not correspond to the wavelet coefficients signature detected by the previously SYM-H 

323 index analysis. The (af J ) 2 (j = 2) presented the highest amplitude coefficients at 03 : 58 UT 

324 and 19 : 52 UT. The maximum amplitude (d-') 2 (j = 3) value occurred at 22 : 40 UT. For (rf 7 ) 2 

325 ( j = 4), the highest wavelet amplitudes occurred at 04 : 16 UT, 05 : 04 UT, 17 : 04 UT, 20 : 00 

326 UT, 21 : 36 UT and 23 : 12 UT. 

327 Figure [5J: shows the wavelet decompositions for the PPT, French Polynesia. As reported 

328 by NOAA, the initial wave arrival could not be determined exactly. For guiding purposes, we 

329 presented the sea level variations at Nuku Hiva, French Polynesia. At all the four decomposition 

330 levels, we detected perceptible (d^) 2 (j = 1,2,3,4) amplitudes between 15 : 00 UT and 20 : 00 

331 UT that might be due to the tsunami propagation. 

332 On 27 December, 2004, the tsunami waves that propagated throughout the Pacific Ocean 

333 arrived at the coast of Peru (see Figure [5};). At Callao, Peru, the initial tsunami wave arrived at 

334 05 : 20 UT with maximum wave height of 0.65 m. The (d-*) 2 for j = 1,2,3,4 decomposition 

335 levels showed what might be tsunami transients between 13 : 00 UT and 19 : 00 UT. 

336 Thus from the methodology applied, the transients in the magnetogram, that are not related 

337 to geomagnetic storms, are identified by the wavelet coefficients. Those transients are not related 

338 to the common geomagnetic disturbance causes, but some remarkable coefficients are related to 

339 the ocean stream. The dashed rectangles in each panel of the figure indicate the identification of 

340 the transients associated to the tsunami magnetic effects. 

341 5.3. 27 -February -20 10 tsunami event 

342 In the region of Callao La Punta, Peru, the observed and computed tsunami time arrival were 

343 coincident, both at 10 : 34 UT with amplitude of up to 0.69 m. We used the sea level measured 

344 at this region as guiding line to the arrival tsunami at HUA and IPM (see Fig|6h-b). Fig|6]is very 

345 similar to Figure|5] also shows the wavelet signatures for the four first decomposition levels ((d J ) 2 

346 where j = 1 , 2, 3, 4). Panels from (a) to (d) correspond to magnetic stations of HUA, IPM, PPT 

347 and TEO, respectively. In each panel, from top to bottom are displayed the corresponding tide 

348 gauge measurement, the SYM-H index, the corresponding magnetogram (Z component), and 

349 the (d-') 2 for j = 1,2,3,4 decomposition levels. On geomagnetic conditions, the 27 February, 

350 2010 corresponded to a very quite day. The Dst index presented the minimum of -2 nT and 

351 the maximum of 4 nT with very smooth variations. The SYM-H also presented very smooth 

352 variations with the minimum of -9 nT and the maximum of 5 nT. 

353 In Fig[6ji, the (d-') 2 for j = 1, 2, 3, 4 decomposition levels presented quite similar behaviors. 

354 Also, they showed a main structure of coefficients which are preceded and followed by a se- 

355 quence of small structures after 1 1 : 00 UT. The highest wavelet coefficient amplitudes occurred 

356 at 1 1 : 26 UT, 15 : 26 UT, between 16 : 42 UT and 18 : 12 UT, 19 : 48 UT and between 21 : 00 

357 UTand21 : 48 UT. 

358 The NOAA Pacific Warning Center predicted the tsunami time arrival at 12 : 05 UT for 

359 the Easter Island, IPM. In Fig |6j), the first decomposition level presented a main structure of 

360 coefficients between 11 : 58 UT and 13 : 24 UT and a secondary pike at 15 : 36 UT, also 
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361 presented a sequence of small structures between 16 : 54 UT and 24 : 00 UT. (d^) 2 (J = 2) 

362 showed the highest coefficients at 12 : 28 UT and 12 : 40 UT and (d j f(j = 3) at 12 : 00 UT. 

363 {d.j) 2 (j = 2,3) presented less structured features than (d-') 2 (j = 1,4). In (d^) 2 (j = 4) was possible 

364 to notice two peaks, one at 12 : 00 UT and the other at 12 : 32 UT, followed by a main structure 

365 of coefficients between 17 : 04 UT and 19 : 28 UT and a secondary structure between 20 : 48 

366 UTand21 : 06 UT. 

367 Fortunately, for PPT, we were able to get the tide gauge measurements at the same location 

368 (see Fig[6p). The observed tsunami initial arrival time was at 17 : 33 UT and the computed time 

369 was at 17 : 47 UT with amplitude up to 0.22 m. The {d^) 2 (j = 1,2, 3, 4) presented surprisingly 

370 similar wavelet signatures. In Fig[6p, the main wavelet coefficient structures are restricted be- 

371 tween 15 : 36 UT and 21 : 52 UT. The tsunami-induced electromagnetic field effects could be 

372 detected about two hours in advance. 

373 We also selected a further located station TEO in Mexico with the purpose of analyzing if 

374 we could detect any geomagnetic tsunami effects. For guiding purposes, we used the tide gauge 

375 measurements at Acapulco. The observed tsunami initial arrival time was at 15 : 47 UT with 

376 amplitude up to 0.66 m. In Fig[6ji, the (d-') 2 (j =1) showed an even more structured signature than 

377 (d J ) 2 (J - 2, 3,4). The highest wavelet amplitudes coefficients for (d^) 2 (j = 1) was at 14 : 26 UT, 

378 16 : 34 UT, 18 : 16 UT, 19 : 40 UT, 20 : 34 UT and 21 : 24 UT. The (d j ) 2 (j = 2, 3, 4) presented 

379 a quite similar signature but less structured features than {d^) 2 (j = 1). It was possible notice the 

380 presence of highest wavelet amplitudes coefficients after 12 : 52 UT for (d^) 2 {j = 2,3,4) that 

381 could be related to the tsunami effects. 



382 5.4. ll-March-2011 tsunami event 

383 On 1 1 March, 201 1, occurred a geomagnetic disturbance which corresponded to a moderate 

384 storm with minimum Dst = -82 nT at 06 : 00 UT and a second energy injection at 18 : 00 UT 

385 with minimum Dst = -67 nT at 22 : 00 UT. So, the abrupt variations of the vertical component 

386 of the geomagnetic field due to the development of this moderate storm will be more emphasized 

387 by the highest amplitudes of the wavelet coefficients. 

388 In Figure|7] the highest coefficients of all the four levels detected in the S YM-H are restricted 

389 to the development of the storm following the second energy injection at 18 : 00 UT, specially 

390 during the main phase of the geomagnetic storm between 18 : 00 UT and 24 : 00 UT presented 

391 the hi ghest wavelet coefficients amplitudes which was usually observed by iDomingues et al 



(2005). All the decomposition levels show quite similar behaviors, showing the highest wavelet 
coefficients (d j ) 2 (j =1,2,3,4) soon after the recovery time of the first storm and associated with 
the fluctuations of the SYM-H index as the second magnetic storm came to the minimum value. 

Figure [8] is very similar to Figure [5] and Fig [6] It also shows the wavelet signatures for the 
four first decomposition levels, (d->) 2 for j - 1,2,3,4. Panels from (a) to (d) correspond to 
magnetic stations of HUA, IPM, PPT and TEO, respectively. In each panel, from top to bottom 
are displayed the corresponding tide gauge measurement, the SYM-H index (both have been 
included for comparison purposes only), the corresponding magnetogram (Z component), and 
the (d-') 2 for j = 1,2,3,4 decomposition levels. We decided to exclude the period between 
18 : 00 UT and 24 : 00 UT, as shown in Figure [7] because the wavelet coefficients in this period 
were related to the abrupt variations caused by the second energy injection in the geomagnetic 
storm. 

In the city of Hanasaki, near to the magnetic station of MMB, the initial phase of the tsunami 
started at 06 : 38 UT and reached the maximum height (0.74 m) at 06 : 57 UT. We used the 
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406 tide gauge measurements at Hanasaki in Figure^ for guiding purposes. The {d^) 2 (j = 1,2,3,4) 

407 showed a very structured signature. Also, (d-') 2 (j - 1, 2, 3) showed a main structure of wavelet 

408 coefficients that could be noticed between 06 : 00 UT and 08 : 00 UT, a secondary structure 

409 between 10 : 00 UT and 1 1 : 20 UT and a smaller structure between 16 : 00 UT and 18 : 00 UT. 

410 The (d^) 2 (j = 4) showed less structured features than {d^) 2 {j = 1, 2, 3) and presented the highest 

411 amplitudes at 07 : 12 UT, 10 : 56 UT and 17 : 04 UT. We can observe at the first decomposition 

412 level, a presence of wavelet coefficients at 06 : 26 UT which corresponds to the tsunami initial 

413 phase. At 06 : 57 UT, when tsunami reached the maximum height, we also notice the presence 

414 of an increase in the coefficient amplitudes with the maximum at 07 : 04 UT higher than the 

415 coefficient at the initial phase. Meanwhile, (d^) 2 (j = 2, 3,4) associated to the tsunami presented 

416 more spread on time and did not present any remarkable amplitudes corresponding to its initial 

417 phase and maximum height. 

418 Near to the KNY station, in the city of Tosashimizu, the tsunami started at 07 : 5 1 UT and the 

419 maximum height (1.28m) was measured at 16 : 58 UT (see Figure[8p). The wavelet signatures of 

420 (d J ) 2 (J = 1 , 2, 3, 4) show very similar structures. The (d^) 2 (j =1,2) showed the highest wavelet 

421 amplitudes around 06 : 04 UT and 17 : 02 UT. And the (d j ) 2 (j = 3,4) showed the highest 

422 wavelet amplitudes around 06 : 08 UT, 06 : 40 UT, 07 : 44 UT, 10 : 40 UT, 16 : 48 UT and 

423 17 : 52 UT. This japanese magnetic station showed surprisingly similar signatures compared to 

424 Figure [8^. We also observe, in third and fourth decomposition levels, the presence of wavelet 

425 coefficients {dj) 2 (j = 3,4) around the tsunami initial phase and maximum height. 

426 In Figure[HJ:, the observed time of the tsunami's arrival at Pago Bay, Guam around 09 : 16 UT. 

427 The tsunami maximum height occurred at 09 : 58 UT and its amplitude was up to 0.54 m. GUA 

428 magnetic station did not show any coefficient {d^) 2 {j = 1,2,3,4) that represented significant 

429 contribution due to the tsunami. At the same time, we observe the presence of small coefficients 

430 (d j ) 2 (j = 1, 2, 3,4) but still noticeable between 10 : 40 UT and 1 1 : 28 UT. 

431 The predicted time of the tsunami's arrival by NOAA at Cairns (Lat. -16.57° and Long. 

432 145.75°), Australia, would be at 15 : 35 UT Unfortunately, at Cape Ferguson, it was not possible 

433 to determine the exactly tsunami time arrival (see Figure [8]i). The tide gauge measurements 

434 registered the tsunami maximum height at 21 : 00 UT with amplitudes up to 0.1 1 m. The CTA 

435 magnetic station did not show any remarkable wavelet coefficient amplitudes {d^) 2 (j = 1 , 2, 3, 4) 

436 due to the tsunami. 

437 In summary, the wavelet coefficient of the first three levels of decomposition showed a good 

438 time localization of the initial phase and the maximum height of the Japanese tsunami at 1 1th of 

439 March, 2011 and they, also, were locally associated with the geomagnetic variations present in 

440 the signal due to the secondary magnetic fields induced by the movement of electrically conduct- 

441 ing sea-water through the geomagnetic field. The red vertical lines in each panel of the figure 

442 indicate the identification of the first transient associated to the tsunami magnetic effects. 

443 Thus from the several events analyzed the wavelet method for identification of the tsunami 

444 magnetic effects could be validated, characterizing an useful tool for this kind of study in order 

445 to deal with very weak but significant amplitude of wavelet coefficients. 

446 6. Conclusions 

447 Dealing with the geomagnetic data analyzed by discrete wavelet transform using Haar wavelet, 

448 this work aimed to implement a methodology for the direct analysis of magnetograms and to ver- 

449 ify the occurrence of tsunami effects on Z-component of the geomagnetic field. The carried out 

450 analysis can be summarized as follows: 
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451 • The wavelet technique is useful to "zoom in" the localized behavior of the geomagnetic 

452 variations induced by the movement of electrically conducting sea-water through the ge- 

453 omagnetic field; i. e., the identification of magnetic transients related to the tsunamis. As 

454 well, from the analysis of the magnetogram data, it was able to localize the initial phase 

455 and maximum height of the tsunamis in some cases. 

456 • The discrete wavelet transform (Daubechies - dbl) provides a very good local characteri- 

457 zation of the signal and may be the most convenient to represent times series with abrupt 

458 variations or steps, i. e., very small and localized variations as the discontinuities in the 

459 Z-component due to tsunamis. 

460 • The discrete wavelet transform is an alternative way to analyze the global influence of the 

461 tsunamis on the geomagnetic field and it could be used as a sophisticated tool to predict 

462 the tsunami arrival, as an example of the Chilean tsunami, the DWT detected transients 

463 about two hours in advance at Papeete. 

464 The results obtained are encouraging even with a few of events. The present study dealt 



465 with only three tsunami events and few magnetic stations. A future study using more events and 

466 more stations should be very good in order to do a more complete analysis of features. The first 

467 interpretation of the results suggests that discrete wavelet transform can be used to characterize 

468 the tsumanis effects on the geomagnetic field. 
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(c) PPT, 26" 1 December, 2004 
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Figure 5: Magnetograms and the wavelet coefficients (d'f for j = 1,2,3,4 (a-d) for LRM, EYR, PPT and HUA, 
respectively. 
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(a) HUA, 21 th February, 2010 



(b) IPM, 27" ! February, 2010 
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(c) PPT, 27"' February, 2010 



(d) TEO, 27"' February, 2010 
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Figure 6: Magnetograms and the wavelet coefficients (d' f for j = 1, 2, 3, 4 (a-d) for HUA, IPM, PPT and TEO, respec- 
tively. 
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11 "'March, 2011 
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Figure 7: Graphics of SYM-H with the corresponding wavelet coefficients (d^) 1 for j = 1, 2, 3, 4. The panel corresponds 
to a japanese tsumani occurrence at 2011. 
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(a) MMB, 11"' March, 2011 



(b) KNY, 11"' March, 2011 
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(c) GUA, 11 "'March, 2011 



(d) CTA, 11 "'March, 2011 
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Figure 8: Magnetograms and the wavelet coefficients (d j ) 2 for j = 1,2,3,4 (a-d) for MMB, KNY, GUA and CTA, 
respectively. 
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